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We analyze the structure of spin-1 Bose-Einstein condensates in the presence of a homogenous 
magnetic field. We classify the homogenous stationary states and study their existence, bifurcations, 
and energy spectra. We reveal that the phase separation can occur in the ground state of polar 
condensates, while the spin components of the ferromagnetic condensates are always miscible and 
no phase separation occurs. Our theoretical model, confirmed by numerical simulations, explains 
that this phenomenon takes place when the energy of the lowest homogenous state is a concave 
function of the magnetization. In particular, we predict that phase separation can be observed in a 
23 Na condensate confined in a highly elongated harmonic trap. Finally, we discuss the phenomena 
of dynamical instability and spin domain formation. 
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I. INTRODUCTION 

The spin degree of freedom of spinor Bose-Einstein 
condensates [JHI, Q leads to a wealth of new phenomena 
not possessed by single-component (spin-frozen) conden- 
sates. New spin-induced dynamics such as spin waves [2|, 
spin-mixing [J] and spin textures 0, [f| have all been pre- 
dicted and observed. These spin-dependent phenomena 
are possible due to the development of optical traps [6[ 
which trap all spin components, rather than just the 
low-magnetic-field seeking spin states of magnetic traps. 
However, the effect of an additional small non-zero mag- 
netic field on the condensate in these optical traps was 
studied even in the seminal theoretical Q and experimen- 
tal [l[ works. In fact the interplay of spin and magnetic 
field has been at the heart of some of the most impressive 
>erimcnts, including the demonstration of 



spin oscillations Q and spin textures 



spinor BEC exr 
spin domains ] 
and vortices [8|. 

A spin-1 spinor BEC in a magnetic field is subjected 
to the well-known Zeeman effect. At low fields the effect 
is dominated by the linear Zeeman effect, which leads 
to a Larmor precession of the spin about the magnetic 
field which is unaffected by spatial inhomogeneities in the 
condensate [9]. At higher magnetic fields the quadratic 
Zeeman effect becomes important, and this can lead 
to much more dramatic effects in the condensate, such 
as coherent pop ulation exchange between spin compo- 
nents @, [HI Ef and the breaking of the single-mode ap- 
proximation (SMA) [l2| ■ The study of the behaviour of a 
spin-1 condensate in the presence of a homogenous mag- 
netic field began with the work of Stenger et al. [l[ where 
the global ground state in ferromagnetic and antiferro- 
magnetic condensates was found to be free of spin do- 
mains. It was later found that the ground state under the 
constraint of fixed magnetization was significantly differ- 
ent to the global ground state and even broke the SMA in 
a harmonic trap [12[ . The SMA continued to be used for a 
homogenous condensate in a homogenous magnetic field, 
leading to the discovery that both phase-matched and 



anti-phase-matched states could exist in both ferromag- 
netic and antiferromagnetic condensates [l3[ . Dynamical 
instability was found to occur in ferromagnetic conden- 
sates in nonzero ma gne tic field, leading to the formation 
of spin domains [Til Il5l| while no spin domain formation 
was predicted to occur [l4| or observed [l6| in antiferro- 
magnetic condensates. It seemed that spin domains in 
antiferromagnetic condensates were only to be found in 
the presence of inhomogenous magnetic fields [l[ . 

Recently, we have shown that, contrary to the com- 
mon belief, antiferromagnetic spin-1 condensates may ex- 
hibit spin domain formation in a homogenous magnetic 
field [17| , provided the condensate is larger than the spin 
healing length. In fact, we found that for a homogenous 
antiferromagnetic BEC with nonzero magnetization all 
states are unstable. The form of the ground state in this 
case was therefore unknown. 

In this work we resolve the ground state phase dia- 
gram of a spin-1 condensate in the absence of a trap- 
ping potential. We show that the translational symme- 
try of a homogenous BEC is spontaneously broken and 
phase separation occurs in magnetized polar condensates 
if the magnetic field is strong enough. An analogous phe- 
nomenon has been predicted and observed previously in 
binary condensates [H, [l9| . To explain the physics be- 
hind phase separation and determine the conditions for it 
to occur, we analyze the excitation spectrum of the inter- 
nal spin degree of freedom of a homogenous condensate. 
In contrast to spatial excitations of the condensate, which 
have the form of sound waves or phonons @, H, [2(| HH , 
we analyze the case when the spin-dependent energy, but 
not the kinetic energy, is increased with respect to the 
ground state [THE?]]- Next, we show that for a range of 
experimental conditions, it is energetically favorable for 
the system to consist of two separate phases composed 
of different stationary states. We demonstrate numeri- 
cally that this phenomenon can be observed in a polar 
condensate trapped in a harmonic optical potential. 

The paper is organized as follows. SectionHTlintroduces 
the theoretical model of a spin-1 condensate in a homoge- 
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nous magnetic field. In Sec. lIIII we study homogenous sta- 
tionary states of the condensate in a magnetic field and 
calculate the energy of each state in terms of the magne- 
tization and the quadratic Zeeman energy shift. We then 
present the internal spin excitation spectra and bifurca- 
tion behavior of the different spin states. In Sec. IIVI we 
analyze the ground state structure and show that phase 
separation can occur in polar condensates and that the 
behavior without a trapping potential can be used to 
predict the ground states in a harmonic trap. Section fVl 
concludes the paper. 



II. MODEL 

We consider dilute spin-1 BEC in a homogenous mag- 
netic field pointing in the z direction. The mean-field 
Hamiltonian of this system is given by, 



H = I dv E ( w v ^* • v ^ + °i n ^ 



J=-,0,+ 



(1) 



where ip-, ipo, ip+ are the wavefunctions of atoms in mag- 
netic sublevels m — — 1,0, +1, M is the atomic mass, 
V(r) is an external potential and n = J2 n j = S I'AjI 2 
is the total atom density. The asymmetric part of the 
Hamiltonian is given by, 



H a = / dri £ 



c 2 



(2) 



\j=-,o,+ 



where Ej is the Zeeman energy shift for state tpj and the 
spin density is, 

F = (F x ,F y ,F z ) = (^F X J, ftFyji, ftPJ) (3) 

where F x y z are the spin matrices [22[ and if = 
("0+ , ipo , ip- ) . The nonlinear coefHcients are given by 
c = 4vrft 2 (2a 2 + a )/3M and c 2 = ^t?(a 2 - a )/3M, 
where as is the s-wave scattering length for colliding 
atoms with total spin S. The total number of atoms 
and the total magnetization 



N = / ndr 



M= J F z dr = J {n + -n^)dv, 



(4) 
(5) 



are conserved quantities. The Zeeman energy shift for 
each of the components, Ej can be calculated using the 
Breit-Rabi formula [23j 

E± = -^Ekfs (l + 4 v / l±a + a 2 ) T giPsB , 
Eo = -~£ H FS (l + 4^1 + « 2 ) , (6) 



where -Ehfs is the hyperfine energy splitting at zero 
magnetic field, a — (gj + 9j)/xsS/£?hfSj where pg is 
the Bohr magneton, gi and gj are the gyromagnetic ra- 
tios of electron and nucleus, and B is the magnetic field 
strength. The linear part of the Zeeman effect gives rise 
to an overall shift of the energy, and so we can remove it 
with the transformation 

H H +{N + M)E+/2 + (N -M)E-/2. (7) 

This transformation is equivalent to the removal of the 
Larmor precession of the spin vector around the z axis 
[TtI l2~il ]. We thus consider only the effects of the 
quadratic Zeeman shift. For sufficiently weak magnetic 
field we can approximate it by SE = (E + +E- — 2Eq)/2 w 
a 2 £"HFs/16, which is always positive. 

The asymmetric part of the Hamiltonian ^ can now 
be rewritten as 

H a = J dr (-5En + j\F\ 2 ^ = J dme{r), (8) 

where the energy per atom e(r) is given by [ici| 

e = - SEpo + ^|f | 2 = -SEpo + °-f (|fj 2 + m 2 ) , 



|fl| 2 =2po(l - Pa) + 2po y/0- - Po) 2 - m 2 cos(#) . (9) 



We express the wavefunctions as ipj 



/npj exp(i8j) 



where the relative densities are pj = rij/n. We also in- 
troduced the relative phase 9 = 9 + + 9- — 29o, spin per 
atom f = F/n, and magnetization per atom m = f z = 
p + — p-. The perpendicular spin component per atom is 
|fJ 2 = / 2 + / 2 . 

The Hamiltonian |T]) gives rise to the Gross-Pitaevskii 
equations describing the mean-field dynamics of the sys- 
tem 

ih ~0f = [£ + c 2 {n± + n - n T )] ip± + c^oV^ i ( 10 ) 
dtp 

ift-gj- = [C- 8E + c 2 (n+ + n-)] ip + 2c 2 ip + tp-ipQ , 

where £ is given by C — — h 2 V 2 /2M + c n + V(r). 

By comparing the kinetic energy with the inter- 
action energy, we can define a characteristic healing 
length £ = 2-nTij \J2Mcqu and spin healing length £ s = 
2-Khj \J2Mc 2 n. These quantities give the length scales of 
spatial variations in the condensate profile induced by the 
spin-independent or spin-dependent interactions, respec- 
tively. Anal ogously , we define magnetic healing length as 
i B = 2Trh/V2MSE. 

In real spinor condensates, the ao and a 2 scattering 
lengths have similar magnitude. The spin-dependent in- 
teraction coefficient c 2 is therefore much smaller than its 
spin-independent counterpart cq. For example, this ra- 
tio is about 1:30 in a 23 Na condensate and 1:220 in a 
87 Rb condensate far from Feshbach resonances [24j. As 
a result, the excitations that change the total density re- 
quire much more energy than those that keep n(r) close 
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to the ground state profile. In our considerations we will 
assume that the amount of energy present in the system 
is not sufficient to excite the high-energy modes, and we 
will treat the total atom density n(r) as a constant. 



III. HOMOGENEOUS STATIONARY STATES 

First, we investigate the homogenous condensate in the 
case of a vanishing potential, V(r) = 0. We look for 
homogenous stationary solutions in the form 

ipj(r,t) = % /7T- e -^ + ^* + ^ , (11) 

where ps = Con/h is a constant. We thus extend the 
studies of [l3| and [lj]]. These solutions are stationary 
in the sense that the number of atoms in each magnetic 
sublevel is constant in time. The relative phase between 
the components may change in time as long as the phase 
matching condition 

H + +/j,- = 2p , (12) 

is fulfilled [13, HH| ■ Because the symmetric part of the 
hamiltonian in ([I]) is constant, the relevant part of the 
hamiltonian is given by Eq. ©. 

The hamiltonian (fTJ) and GP equations (fTTTj) are invari- 
ant under the gauge transformation ipj — > ipj exp(— iff) 
and rotation around the z axis ipj — > ipj exp(— iF z ^y), 
which transform the wavefunction components according 
to 

/ i>+ \ ( e" i7 ^+ \ 

I _> e -*> \ iPo . (13) 

V i>- ) \ e 1 ^- / 

Hence the solutions can be classified using the rela- 
tive densities pj = nj/n and a single relative phase 
8 = 8 + +8- — 28q, with the chemical potentials pj given as 
solutions to Eqs. (fit))) . We note that for stationary solu- 
tions the relative phase must take one of two values, 8 = 
or 8 = it. We call the former "phase-matched" states 
(PM) and the latter "anti-phase-matched" (APM) states. 
The names derive from the fact that within the contin- 
uum of states satisfying the spin rotations (fTB")) there is 
a set (ip+,ipo,ip-) with all components in phase for the 
PM states, and a set with ip + and ipo in phase but ir out 
of phase with tp_ for the APM states. 

The following analysis is also applicable to nonhomoge- 
neous condensates within the single-mode approximation 
(SMA), which assumes that the spin components share 
the same spatial profile @, 0; [H| , after replacing n with 
(n). This assumption is true eg. when the condensate 
size is much smaller than the spin healing length £ s . 

Stationary solutions of the system ([TO]) may have one, 
two or three nonzero components. We examine each case 
separately and then examine the existence regions and 
bifurcation behavior of the three states together. 



A. One-component solutions (po, P±) 

If only one component has nonzero atom density, we 
have two qualitatively distinct possibilities: 

1. po = 1. This state exists for m = 0. From Eq. © 
and substitution of (| 1 1 [) into (fl"Uj) we find the chem- 
ical potential of the m = component, the perpen- 
dicular spin per atom and the energy per atom, 

hp = -5E, |f±| 2 = 0, e P0 =-SE. (14) 

2. p + = 1 or p_ = 1. These two states exist for 
m = lorrn = — 1, respectively. Following a similar 
procedure to the case above we find, 

HfJ.+ = c or hn-=c, |fi| 2 =0, e p± = - . (15) 

where we have introduced a shorthand notation for 
the effective interaction coefficient, c = c^n. 

B. Two-component solutions (2C) 

Here one can in general choose the vanishing compo- 
nent arbitrarily, but only one choice turns out to be a 
stationary solution. 

1. po = 0. One distinct state exists for any values of 
m, c, and SE. 

hp,±=±cm, |fjj 2 = 0, e 2 c = ^m 2 . (16) 

2. p + = or p_ = 0. Due to the spin-dependent 
interaction these cases are non-stationary leading 
to generation of the third component, as is evident 
from the final terms in Eqs. (fTU]) . 

C. Three-component solutions (PM, APM) 

We can derive the relationship between 5E/c 7 m and 
po from the phase matching condition (fl"2"|) and the GP 
equations (fTU)) . 

SE (l-p )(l-2p )-m 2 

— = 1 - 2p + s == — , (17) 

c ^(l - po) 2 - m 

where s = 1 for PM states and s = — 1 for APM states. 
This condition can be alternatively derived from the en- 
ergy functional as de/d{po, d)\ m ,SE,n = 0, since the sta- 
tionary states correspond to energy extrema under con- 
straint on m and n. Note that SE/c can be positive or 
negative depending on the sign of c 2 . 

It can be shown that if m ^ 0, there can be at most 
one distinct PM and one APM solution from the interval 
Po £ [0, 1] for a given SE/c and m. More specifically we 
find the following: 
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1. PM solutions (0 = 0) 

(a) For ra/0, there is one PM solution for each 
SE e (— oo, 6Epm), where we find <5£pm = 
c (1 + Vl - m 2 ) from |T7J) for p = 0. 

(b) For m = 0, the solutions exist in the interval 
SE e (-2c, 2c). 

2. APM solutions (9 = n) 

(a) For m 7^ 0, there is one APM solution for 
each SE £ (SEapm,+oo), where <5apm = 
c (l — Vl — wi 2 ) (again found from (fTT|) with 
Po - 0). 

(b) For m = 0, an infinite number of degener- 
ate solutions exist at SE = (no magnetic 
field), for any value of p n e [0,1], however 
APM states do not exist at nonzero field. 

We see from ([9]) that for three-component solutions f± 
is nonzero for both PM and APM states, hence these 
solutions break the U(l) rotational symmetry. In fact, 
the state investigated in [2l[ is identical to the PM state, 
which is the ground state of a ferromagnetic condensate 
whenever m ^ or SE < 2\c\. Moreover, as the chemi- 
cal potentials of the three components arc in general not 
equal the spin vector rotates around the z axis, in addi- 
tion to the Larmor rotation, at a rate proportional to the 
difference in chemical potentials. 

To obtain further information about the nature of the 
three-component solutions we examine some of the lim- 
iting behavior. By expanding p$ about 1 — \m\ we can 
determine the behavior in the limit \SE/c\ — *■ +oo, 



p = (l- \m\) - 



2SE- 



M(i-H) 2 , 



/i = 2|m|(l-|m|), 



(18) 



c\m\ 1 — 



3(A)pm = -SE(1 - \m\ 

hp = -SE + 0(c) , hp± = -SE[1 T sgn(m)] + 0(c) 



|m| 
2 



On the other hand, when SE — > SE^ PM (close to the 
bifurcation point), 



Po 

fl 
hpo 



1 \5E-SE {A)PM \ 



SE, 



(A)PM 



(19) 
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FIG. 1: Diagram of existence of two- (2C) and three- 
component (PM or APM) homogenous stationary states in 
spin-1 condensates. In addition to the solutions shown on 
the diagram, one-component solutions pj = 1 exist with 
j = — ,0, +. The dotted and dashed lines at M = indi- 
cate the absence of a PM or APM state respectively. 



of the diagram corresponds to polar condensates (c2 > 0) , 
and the lower half describes ferromagnetic BECs (c2 < 
0). There is clearly a region of coexistence of PM and 
APM solutions in polar condensates. For ferromagnetic 
condensates APM states only exist at zero magnetic field 
(and zero magnetization). 

The order parameter space of the two- and three- 
component solutions in nonzero magnetic field is U(l) x 
U(l) (toroid), due to the two symmetries, gauge and 
rotation around z axis, which leave the phase 9 and 
the atom density in each of the components unchanged. 
However, the one-component solutions have only the 
U(l) parameter space (the same as for a scalar conden- 
sate), due to the equivalence of gauge transformation and 
rotation. When no magnetic field is present, the situation 
is different because PM and APM states become degen- 
erate for all values of po and form the families of polar 
and ferromagnetic states together with p± and po states, 
respectively (see Fig. [4j . The order parameter manifolds 
of these families are SO(3) for the ferromagnetic [2| and 
U(l) x S 2 /Z 2 for the polar states family [23, HI]. 



E. Internal spin excitations spectra and 
bifurcations 



D. Existence and order parameter spaces 

Combining the above results, the complete diagram 
of existence of the two-component and three-component 
states is shown in Fig.Q] Because SE > 0, the upper half 



In scalar condensates, where only one spin component 
is present, the Bogolubov theory can be used to describe 
spatial excitations of the condensate, which have the form 
of sound waves or phonons [20]. In spinor condensates, 
modes of a similar nature have been studied [1, H, Hl|] , but 
another degree of freedom is also available. One can con- 
sider internal spin excitations, where the spin-dependent 
energy, but not the kinetic energy, is increased with re- 
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(a) Ferromagnetic condensate 

-l> / |m| 

5E/c -0.5 \/ 

.(. 00 

(b) Antiferromagnetic condensate 
PM — 2C 




FIG. 2: Normalized energy per atom e/(\c2\n) as a function 
of magnetization |m| = \A4\/N and normalized quadratic 
Zeeman energy SE/(c2n) for homogenous stationary states 
of spin- 1 condensates, including two-component (2C), phase- 
matched (PM) and anti-phase-matched (APM) states. The 
single-component p — 1 state is also shown, (a) Ferromag- 
netic condensate; (b) Anti-ferromagnetic condensate, with 
lower panel showing a zoom of the upper panel. 



spect to the ground state. In contrast to the phonon-type 
excitations, the excitation spectrum of the internal spin 
states is discrete. The stationary solutions, described in 
the previous section, form a set of such modes when a 
magnetic field is applied to a homogeneous condensate. 
We argue that the energy of the system can be exchanged 
between the spin modes and spatial excitations, which 
has important consequences for the condensate dynam- 
ics. 
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FIG. 3: (a,b) Normalized energy per atom e/(|c2|n) in func- 
tion of normalized quadratic Zeeman energy SE/(c2n) for ho- 
mogenous stationary states in the M = case for (a) ferro- 
magnetic and (b) polar condensate. The arrows in (a) show 
the scenario of the quenched condensate experiment [§]. (c,d) 
Energy in function of the absolute value of the magnetiza- 
tion |m| for (c) SE/(c 2 n) = 0.5 and (d) 5E/(c 2 n) = 1.5. The 
black dots indicate bifurcation points where SE = <5-Eapm and 
SE — SEpm, respectively. The energy of the APM state is a 
slightly concave function of the magnetization, which gives 
rise to the phase separation (see Sec. IIV[) . 



The dependence of the energy per atom given by 
Eq. @ on magnetic field and magnetization for the spin 
states studied in the previous section is shown on three- 
dimensional plots in Fig.O Note that renormalized vari- 
ables e/(\c2\n) and bEj{cin) are used, which allows us 
to include all possible configurations of spin-1 condensate 
just in two universal graphs (there is no fixed parameter). 

The energy dependence cross-section in the particular 
case m = is shown in Fig. (2Ua,b). It is clear from 
Fig. that when the Zeeman energy is decreased, the 
po state ceases to be the lowest-energy state of the fer- 
romagnetic condensate, and the PM state becomes the 
ground state [2lJ. This fact has been utilized in the ex- 
periment by Sadler et. al. 0], where the ferromagnetic 
BEC prepared in the po state in strong magnetic field was 
suddenly quenched to the low magnetic field regime. As 
the condensate relaxed locally to the PM state in a con- 
servative process, an amount of energy was transformed 
from the spin to the kinetic energy, which allowed the for- 
mation of spin domains and topological excitations. As 
we will show later, this excess energy was necessary since 
no spin domains exist in the ground state of a ferromag- 
netic condensate. In Fig. [3Jc,d) we show dependence of 
energy on magnetization for fixed magnetic field strength 
and indicate the points of bifurcation. 

We summarise the degeneracies and bifurcations be- 
tween the various spin states in the following list: 

1. The three-component (PM, APM) and two- 
component (2C) solutions become identical to one- 
component solutions p± at m = ±1. 
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FIG. 4: Schematic picture of equivalences and bifurcations of 
homogenous stationary states in various regimes of 8E/(c2n). 
Solid lines indicate that two states become identical for a 
particular value of |mj, and dashed lines correspond to bi- 
furcations occuring at a given value of 5E. In the SE — 
case (no magnetic field), the two circles show that degen- 
erate states form the families of polar and ferromagnetic 
states, which are the ground states of polar and ferromag- 
netic condensates, respectively @|. Note that for ferromag- 
netic condensates(c2 < 0), the energy grows to the left. 



2. The PM states bifurcate from two-component so- 
lutions (2C) at SE = SEpm- The APM states bi- 
furcate from 2C solutions at 8E — SEapm- 

3. The APM states become identical to po states at 
to = 0. 

4. In a ferromagnetic condensate (c < 0), PM states 
become identical to the po state at m = if 
SE > 2\c\. However, for SE E (0,2|c|) there ex- 
ists a separate PM state with to = which is not 
equivalent to po, see Fig. [3Ja). Hence SE = —2c, 
to = is a bifurcation point. 

5. At SE = (no magnetic field), all the APM states 
have |f| = 777 = and are degenerate and con- 
tinuously connected to the po state. These states 
together form the family of polar states, which is 
the ground state of polar condensates 0, [27| • 

6. At S E = 0, all the PM states (with different values 
of to) have |f | = 1 and become degenerate and con- 
tinuously connected to the p± states. These states 
form the family of ferromagnetic states, which are 
the ground states of ferromagnetic condensates [1] . 

These connections between states, together with the 
energy hierarchy, are schematically collected in Fig. 2] 
Equivalence between two states at to — or \m\ = 1 is 



indicated with continuous lines, while bifurcations (oc- 
curring with either changing to or SE/(c2n)) are marked 
by dashed lines. For example, 2C states become equiva- 
lent to p± states (i.e. either p+ or p_) whenever m — > ±1. 
On the other hand, APM states bifurcate from 2C states 
when the quadratic Zeeman energy crosses the value 
SEapm, and separate APM and 2C states exist for the 
same value of to for SE > SEapm- The two circles in the 
middle correspond to the polar and ferromagnetic state 
families at zero magnetic field. 



IV. GROUND STATES AND PHASE 
SEPARATION 

A. No trapping potential 

The ground states of spin-1 condensates in homoge- 
nous magnetic field have been studied in a number of 
previous works [l], [H, HH, [H, [2?| • The most common pro- 
cedure [l|, l2l| involves minimization of the energy func- 
tional with constraints on the number of atoms N and the 
total magnetization M . The resulting Lagrange multipli- 
ers p and q serve as parameters related to the quadratic 
Zeeman shift SE and the magnetization to. An alterna- 
tive method, elaborated in pjj, consists of minimization 
of the energy functional in the parameter space of phys- 
ically relevant variables B and to. Most of the previ- 
ous studies, however, were assuming that the condensate 
remains homogenous and well described by the single- 
mode approximation; in particular, the spatial structure 
observed in [l[ resulted from the applied magnetic field 
gradient, but the BEC was assumed to be well described 
by the homogenous model at each point in space. In 
Ref. [l2j , the breakdown of the single-mode approxima- 
tion was shown numerically for a condensate confined in 
a harmonic potential. 

We correct the previous studies by showing that when 
the condensate size is larger than the spin healing length 
£ s , the translational symmetry is spontaneously broken 
and phase separation occurs in magnetized polar conden- 
sates if the magnetic field is strong enough. This phe- 
nomenon takes place when the energy of the spin state 
with the lowest energy is a concave function of to for 
a given SE. On the contrary, the energy is always a 
convex function of to for the ferromagnetic condensate, 
and no phase separation occurs. Note that phase separa- 
tion has been previously predicted in binary condensates 
[l8[ and in ferromagnetic condensates at finite tem- 
perature [29j . 

Two types of domain structures, depicted in Fig. [5l 
are composed of two different stationary states connected 
with a shaded region where all three components are 
nonzero to ensure proper matching of the chemical poten- 
tials fTS]) . These two domain states have the advantage 
that the perpendicular spin is nonzero only in the transi- 
tory region, hence their energy is relatively low in polar 
condensates. In fact, these are the only phase separated 
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FIG. 5: Schematic structure of the phase separated states 
a) p± + po an b) 2C + po. The shaded region, in which all 
three components are nonzero, has the approximate extent 
of one spin healing length £ s or magnetic healing length £g, 
whichever is greater. The relative size of the domains is in- 
dicated with arrows, (c, d) The corresponding wavefunction 
profiles obtained numerically with periodic boundary condi- 
tions in the case of 23 Na for m = 0.5 with (c) SE/ (c2n) = 0.8 
and (d) 8E/(c2n) — 0.23. The n+, no, and n_ components 
are depicted by dash-dotted, dashed, and dotted lines, respec- 
tively. The solid lines show the total density. 



Condensate 


Parameter range 


Ground state 


Ferromagnetic 


2 < -M- and m = 

— \C2\n 


PO 




j^- < 2 or m / 


PM 


Polar 


m = 


PO 




SE <- m 2 
C2TI — 2 


2C 




^<^<|andm^0 


2C + po 




1 < iS. and m / 


p± + po 



TABLE I: Ground states of spin-1 condensates in homogenous 
magnetic field. The states 2C + po and p± + po correspond to 
phase separation (see Fig. [5} . 



states that can be the ground states of a homogenous 
condensate. Their energies per atom in the limit of in- 
finite condensate size, which allows for neglecting of the 
relatively small intermediate region are 



2p±+Po = H e P± + C 1 - H)e P0 , 



e 2C+p 



m 2C 







+ 


(i- m ) 


1 m—m2C, 


V rn 2C J 



(20) 



where rn = M/N is the average magnetization and the 
magnetization of the 2C component rri2C is a free param- 
eter that has to be optimized to obtain the lowest energy 
state. 

The ground states can be determined by comparing 
energies of the phase separated states with the energies 
of the homogenous solutions of Sec. IHII The results for 
both polar and ferromagnetic condensates are collected 
in Table HI In the cases when no phase separation occurs, 
our results are in agreement with those obtained in 
Note that we assumed that the condensate size is much 
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FIG. 6: Ground state phase diagram of the polar condensate. 
The symbols correspond to numerical data obtained for the 
parameters of 23 Na, with solid triangles representing 2C, open 
circles 2C+po and open squares p± + po- The solid lines and 
shading are given by the analytical formulas from Table [I] 



larger than £ s and £b. For small condensates, the results 
of [12| are correct. 

In the case of high magnetic field strength, one of 
the Zeeman sublevels is practically depleted [12| and 
the condensate becomes effectively two-component. The 
existence of the p± + po phase in a polar condensate 
can then be understood within the binary condensate 
model fl8l. ITfjj . We note that the experiment reported in 
Ref. [301 ] . performed in this regime, can be viewed as the 
first confirmation of phase separation in spin-1 BEC in a 
homogenous magnetic field. However, the ground state 
was not achieved, and a multiple domain structure was 
observed. 

In Fig. [6] we present the phase diagram of polar con- 
densates, obtained both numerically and using analytical 
formulas from Table |U The ground state profiles for a 
quasi- ID condensate were found numerically by solving 
the ID version of Eqs. E3 



ih 



ih 



dip 



± 



at 

dtpo 
dt 



C + c 2 (n± + h - n T ) ip± + c 2 ip ip^ , (21) 
C — SE + c 2 (h + + ri_) ip + 2c 2 tp + ip-ipo , 



with C = —(h 2 /2m)d 2 /dx 2 + en, where en = 47kjj_(2a 2 + 
ao)/3, c 2 = 4hio±(a 2 — ao)/3, J dx^ = N and lu± is 
the transverse trapping frequency. We imposed periodic 
boundary conditions on ipj(x) and used the parameters 
corresponding to a 23 Na BEC containing N — 5.2 x 10 4 
atoms confined in a transverse trap with frequency lo±_ = 
2tt x 10 3 . The Fermi radius of the transverse trapping 
potential is smaller than the spin healing length, and the 
nonlinear energy scale is much smaller than the trans- 
verse trap energy scale, which allows us to reduce the 
problem to one spatial dimension [24], [3l|. The solu- 
tions were found numerically using the normalized gradi- 
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x [|im] 



FIG. 7: Ground state profiles in a harmonic trap potential. 
Phase separation occurs in the polar 23 Na condensate when 
the magnetic field strength is increased from (a) B = 0.1G, 
5£/(c 2 n max ) = 0.09 to (b) B = 0.12G, 5E/(c 2 n max ) = 0.13 
and (c) B — 0.25G, 5-E/(c2n m ax) = 0.56. For comparison, 
the ground state of a 87 Rb condensate is shown in (d) for 
B = 0.2G, 5E/(c*n max ) = -0.41. The n+, no, and n_ 
components are depicted by dash-dotted, dashed, and dot- 
ted lines, respectively. The solid lines shows the total density. 
Other parameters are N — 2.1 x I0 3 , = 2n x 70 ( 23 Na), 
W|| = 2tt x 48 ( 87 Rb), cjx = 2n x I0 3 and m = 0.5. 



ent flow method [32j j . which is able to find a state which 
minimizes the total energy for given N and M, and ful- 
fills the phase matching condition (fT2")l . The stability of 
the resulting states was verified using numerical time evo- 
lution according to Eqs. (f2l"j) . The slight discrepancy be- 
tween numerical and analytical results can be accounted 
for by the finite size of the condensate (the box size was 
~ 10 £ s ), and by the deviation from the assumption that 
the total density is constant (see the discussion at the end 
of Sec. HI| . Due to the finite value of the ratio cz/cq there 
is a slight density modulation, as is evident in Fig.[5^c,d). 



B. Condensate trapped in a harmonic optical 
potential 

The results from the preceding subsection can be ver- 
ified experimentally in configurations involving toroidal 
or square-shaped optical traps (33|. However, in most ex- 
periments on BECs, harmonic potentials are used. The 
relevance of these results is not obvious in the case of 
harmonic trapping, since the coefficient SE/fan), one of 
the main parameters controlling the condensate proper- 
ties, varies in space due to the varying total density n. 

The ground states in a highly elongated harmonic trap, 
where the parallel part of the potential has the form 
V(x) = iMwjja; 2 , are presented in Fig. [7] We can see 
that as the magnetic field strength is increased, phase 
separation occurs and the p±+po domain state is formed. 
However, in contrast to the previous case, the transi- 
tion is not sharp, and in particular there is no distinct 
2C + pq phase for any value of the magnetic field. Note 



that the state in Fig. [TJa) is also spatially separated due 
to different Thomas-Fermi radii of the and ip + com- 
ponents; however, this is an example of potential separa- 
tion, as opposed to phase separation [l9| , since it is does 
not occur in the absence of the potential. On the other 
hand, Fig. [7]Jd) shows that the components of ferromag- 
netic condensate are miscible even in the regime of strong 
magnetic field. In the regions where the wavefunctions 
overlap, the relative phase is equal to 9 = for ferro- 
magnetic and 9 = tt for polar ground states, since these 
configurations minimize the spin energy 

The characteristic feature of phase separation in the 
polar BEC is that the m = domain tends to be local- 
ized in the center of the trap, as shown in Fig. EJb) and 
(c) . This can be explained by calculating the total asym- 
metric energy of the condensate ([8} , again assuming that 
the contribution from the intermediate region connecting 
the domains is negligible, 



ff,w / drn(-SE)+ f drn^ 

Jpo Jp ± 

= -SE(N-\M\) + ^(n) p± , 



(22) 



where 



is the mean condensate density within the 



area of the p± domain. We see that the energy will be 
the lowest if this domain is localized in the outer regions, 
where the condensate density is low. 



C. Spin domain formation 

Our results presented above show that the domain 
structure forming in polar condensates is absent in ferro- 
magnetic BECs. This may seem to contradict the com- 
mon understanding of ferromagnetism and the results of 
the quenched BEC experiment Ref. [8| . The conventional 
picture of a ferromagnet involves many domains pointing 
in various directions separated by domain walls. Similar 
structure has been observed in Ref. &]. However, these 
cases correspond to the situations when there is an ex- 
cess kinetic energy present in the system, due to finite 
temperature or excitation of the spatial modes. On the 
other hand, our study is limited to the ground states at 
T = 0. It is easy to see from Eq. © that in zero magnetic 
field the ground state of a ferromagnetic BEC will always 
consist of a single domain with maximum possible value 
of the spin vector |f | = 1, pointing in the same direction 
at all points in space. However, when the temperature is 
finite, more domains can be formed each with a different 
direction of the spin vector. 

We emphasize that the domain structure of the ground 
state in polar condensates is very different from the do- 
mains formed when the kinetic energy is injected in the 
system as in Ref. 0. The latter constantly appear and 
disappear in a random sequence d, 0, [H, [l?], [H, Hj| . 
On the contrary, the ground state domains are stationary 
and are positioned in the center of the trap. They exist 
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in the lowest-energy state, while the dynamical domains 
require an amount of kinetic energy to be formed. The 
ground state domains can be prepared in an adiabatic 
process, involving adiabatic rf sweep or a slow change 
of the magnetic field [30lj33|. while the kinetic domains 
require a sudden quench [a, l34j . 



the magnetic field strength. The time required for the 
development of instability may be much longer than the 
condensate lifetime 0. 



V. CONCLUSIONS 



D. Dynamical stability 

The dynamical instability of ferromagnetic conden- 
sates that leads to spontaneous formation of s pin do- 
mains has been investigated theoretically [3, HU, HH and 
observed in experiment |8j. An analogous phenomenon 
has been predicted recently for polar condensates in pres- 
ence of magnetic field [l7|. Here we correct the results 
of Ref . [ItJ , by noting that the po = 1 state is stable in 
ferromagnetic condensates for SE > 2\c2\n, and the 2C 
(po = 0) state is stable in polar BECs if SE < m 2 /2. 
Both states become the ground states for these values of 
parameters. By investigating stability in various ranges 
of parameters, we are able to formulate a phenomenologi- 
cal law governing the dynamical stability of condensates: 

1. The only stable state for both polar and ferromag- 
netic BECs in finite magnetic field is the ground 
state, as shown in Table fl] 

2. In zero magnetic field, the same is true for fer- 
romagnetic condensates; However, all stationary 
states of polar condensates are dynamically stable 
in zero magnetic field [3, [l?], [36| . 

The reason for the stability of polar condensates in van- 
ishing magnetic field case is not yet clear. We note that 
the polar condensates in weak magnetic field may also 
be effectively stable on a finite time scale. As shown in 
Ref. [13], in this latter case the instability growth rate 
of unstable modes is proportional to the fourth power of 



We have studied the ground state of a spin-1 BEC 
in the presence of a homogenous magnetic field with 
and without an external trapping potential. We have 
found that without a trapping potential the translational 
symmetry can be spontaneously broken in polar BEC, 
with the formation of spin domains in the ground state. 
We have determined the ground-state phase diagram in 
the space of magnetization versus magnetic field divided 
by density, and demonstrated the different phases, each 
characterized by the type of nonvanishing components. 
We have found good agreement between the numerical 
calculation of the phase diagram and the analytical pre- 
dictions based on the homogenous states. We have shown 
that these results may be used to understand the ground 
state structure in the presence of a trapping potential by 
mapping the locally varying density in the trap to the 
homogenous state. We have found that, depending on 
the magnetic field, the antiferromagnetic BEC ground 
state in the trap displays pronounced spin domains for 
a range of possible experimental conditions. Finally, we 
have discussed the relationship between the phenomenon 
of phase separation and the dynamical instability leading 
to the formation of dynamic spin textures. 
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